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Foreword. This fundamental memoir contains the ORIGINAL form of the solution of Schwarzschild's problem. 
It is regular in the whole space-time, with the only exception of the origin of the spatial co-ordinates; consequently, 
it leaves no room for the science fiction of the black holes. (In the centuries of the decline of the Roman Empire 
people said: "Graecum est, non legitur" ...). 

§1. In his work on the motion of the perihelion of Mercury (see Sitzungsberichte of November 
18th, 1915) Mr. Einstein has posed the following problem: 
Let a point move according to the prescription: 



S J ds = 0, 

where ^ 



. ds = yjY^g^dx^dxv fi,u = l, 2, 3, 4, 

where the g^ v stand for functions of the variables x, and in the variation the variables x must be 
kept fixed at the beginning and at the end of the path of integration. In short, the point shall move 
along a geodesic line in the manifold characterised by the line element ds. 

The execution of the variation yields the equations of motion of the point: 



d 2 x a _ s-^ dx^ dx v a g-123 4 (2) 

ds 2 ~ 2^ V ds ds > 1,2,3,4, w 

where 

V ~ lY \dx v + dx, dxj> (3) 

and the g a @ stand for the normalised minors associated to g a p in the determinant \g^ v \. 

According to Einstein's theory, this is the motion of a massless point in the gravitational field 
of a mass at the point X\ = X2 = x% = 0, if the "components of the gravitational field" V fulfil 
everywhere, with the exception of the point X\ = X2 = x% = 0, the "field equations" 
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and if also the "equation of the determinant' 



1 9v, 



= -1 



(5) 



is satisfied. 



The field equations together with the equation of the determinant have the fundamental prop- 
erty that they preserve their form under the substitution of other arbitrary variables in lieu of xi, 
X2, xs, X4, as long as the determinant of the substitution is equal to 1. 

Let Xi, X2, X3 stand for rectangular co-ordinates, x 4 for the time; furthermore, the mass at 
the origin shall not change with time, and the motion at infinity shall be rectilinear and uniform. 
Then, according to Mr. Einstein's list, loc. cit. p. 833, the following conditions must be fulfilled 
too: 

1. All the components are independent of the time £4. 

2. The equations g p 4 = g± p = hold exactly for p = 1, 2, 3. 

3. The solution is spatially symmetric with respect to the origin of the co-ordinate system in the 
sense that one finds again the same solution when x\, X2, X3 are subjected to an orthogonal 
transformation (rotation). 

4. The vanish at infinity, with the exception of the following four limits different from zero: 



The problem is to find out a line element with coefficients such that the field equations, the equation 
of the determinant and these four requirements are satisfied. 

§2. Mr. Einstein showed that this problem, in first approximation, leads to Newton's law 
and that the second approximation correctly reproduces the known anomaly in the motion of the 
perihelion of Mercury. The following calculation yields the exact solution of the problem. It is 
always pleasant to avail of exact solutions of simple form. More importantly, the calculation proves 
also the uniqueness of the solution, about which Mr. Einstein's treatment still left doubt, and 
which could have been proved only with great difficulty, in the way shown below, through such an 
approximation method. The following lines therefore let Mr. Einstein's result shine with increased 
clearness. 

§3. If one calls t the time, x, y, z, the rectangular co-ordinates, the most general line element 
that satisfies the conditions 1-3 is clearly the following: 



where F, G, H are functions of r = \J x 2 + y 2 + z 2 . 

The condition (4) requires: for r = 00 : F = G = 1, H = 0. 

When one goes over to polar co-ordinates according to x = r sin # cos <f>, y = r sin -& sin cfi, z = 
rcos-tf, the same line element reads: 



Now the volume element in polar co-ordinates is equal to r 2 sin 'ddrdddcj), the functional determinant 
r 2 sin ■d of the old with respect to the new coordinates is different from 1 ; then the field equations 



544 = 1, 5n = 522 = 533 = -1- 



ds 2 = 



Fdt 2 - G{dx 2 + dy 2 + dz 2 ) - H{xdx + ydy + zdz) 2 



ds 2 = Fdt 2 - G(dr 2 + rW + r 2 sin 2 M$ 2 ) - Hr 2 dr 2 
= Fdt 2 - (G + Hr 2 )dr 2 - Gr 2 {di) 2 + sin 2 M<t> 2 ). 



(6) 
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would not remain in unaltered form if one would calculate with these polar co-ordinates, and one 
would have to perform a cumbersome transformation. However there is an easy trick to circumvent 
this difficulty. One puts: 



Xi 



3 



X2 = — COS $, Xs 



(7) 



Then we have for the volume element: r 2 dr sin Ddfid^ = dx\dx-2.dx%. The new variables are then 
polar co-ordinates with the determinant 1. They have the evident advantages of polar co-ordinates 
for the treatment of the problem, and at the same time, when one includes also t = x±, the field 
equations and the determinant equation remain in unaltered form. 
In the new polar co-ordinates the line element reads: 



ds 2 



Fdxl 



- d +^)dx 2 -Gr 2 



dx\ 2 



\-x\ 



+ dxi(l-x£ 



(8) 



for which we write: 



ds 2 = f 4 dx 2 - hdx{ - h—^ - hdx%(\ - x 2 ). (9) 

1 x 2 

Then /i, f 2 = fe, f 4 are three functions of X\ which have to fulfil the following conditions: 

1. For xi = oo : f 1 = 1/r 4 = (3xi)- 4 / 3 , f 2 = f 3 = r 2 = (3xi) 2 / 3 , / 4 = 1. 

2. The equation of the determinant: /i ■ f% ■ fs ■ = 1. 

3. The field equations. 

4. Continuity of the /, except for x\ = 0. 

§4. In order to formulate the field equations one must first form the components of the 
gravitational field corresponding to the line element (9). This happens in the simplest way when 
one builds the differential equations of the geodesic line by direct execution of the variation, and 
reads out the components from these. The differential equations of the geodesic line for the line 
element (9) result from the variation immediately in the form: 



f d 2 Xl idfWdxA 2 idf 1 fdx 1 \ 2 ldfc 

Jl ds 2 2 d Xl \ds J 2 8 Xl \ds J 2 dx 1 



- ^ X2 + ®f 2 1 dxi dx2 + /2^2 ( dx 2 \ 2 + dx^ 2 

1 — x\ ds 2 dxi 1 — x\ ds ds (1 - x|) 2 \ ds J \ ds 



2 ,d 2 x 3 df 2 



2 dxi dx 3 



= / 2 (l-x^ + ^(l-^-- 



dx2 dx 3 
ds ds 



q j. d 2 x 4 ^ 9/4 dx\ dx4 
ds 2 dx\ ds ds 

The comparison with (2) gives the components of the gravitational field: 
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11 " 2/xas!' 122 + 2f 1 dx 1 l-x* 2 
i33 " + 2/ 1 0x 1 (1 X2h 

1 44 



r 2 - X2 r 2 -- x m- x 2 ) 

1 21 — o -f a™ ' 1 22 — 1 2' 33 — X 2(,J- X 2 j, 



2f 2 dx 



i31 " 2^d*7' i 32-+T— 

1 1 9/4 



1 41 



2 / 4 ax a 



(the other ones are zero). 
Due to the rotational symmetry around the origin it is sufficient to write the field equations only 
for the equator [x2 = 0); therefore, since they will be differentiated only once, in the previous 
expressions it is possible to set everywhere since the beginning 1 — x\ equal 1. The calculation of 
the field equations then gives 

d(\dj x \ 1/19/jV /1%\ 2 + 1/1^ 2 



dx 1 \f 1 dx 1 J 2\f 1 dx 1 J \hdx 1 J 2\f 4 dx 

„) JL fi»M_ 2+ 1 



dx 1 \f 1 dx 1 J " hh\dxi 
d ( 1 dh\ 1 (dh 
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dx 1 V h 9xi J /1/4 V dx\ , 
Besides these three equations the functions f\, f 2 , f% must fulfil also the equation of the determinant 

, Wf2f , . 1 dh 2 df 2 1 a/ 4 

d) /1/2/4 = 1, *■ e. -rj^— + + rr7T~ = °- 

h 0x1 h oxi h dxi 

For now I neglect (b) and determine the three functions /1, / 2 , / 4 from (a), (c), and (d). (c) can 
be transposed into the form 

c „ 8(1 df 4 \ _ 1 dh dh 



dxi V h dxi J hf 4 dxi dxi 
This can be directly integrated and gives 

1 d f 

c") = — — — = afi, (a integration constant) 
The addition of (a) and (c') gives 



9xi\/i3xi hdxxj \hdxxj 2\f 1 dx 1 f 4 dx 
By taking (d) into account it follows 



-2 



d ( 1 df 2 \ J 1 df 2 ^ 1 



dx 1 \f 2 dx 1 J \f2dx1 
By integrating 

1 3 p 

1 ^ = -xi + — (p integration constant) 

h dx i 

or 

1 9/ 2 2 



f 2 dx 1 3x 1 +p' 

By integrating once more, 

f 2 = A(3xi + p) 2 ^ 3 . (A integration constant) 
The condition at infinity requires: A = 1. Then 

/ 2 = (3x 1+ p) 2 / 3 . (10) 

Hence it results further from (c") and (d) 

dfi a a 

= a/1/4 - 



By integrating while taking into account the condition at infinity 

/ 4 = l-a(3x 1 + / 9)- 1 / 3 . (11) 

Hence from (d) 

(3x 1+ p)- 4 / 3 

/1_ l-a(3x 1+/ o)-V3- W 

As can be easily verified, the equation (b) is automatically fulfilled by the expressions that we found 
for /1 and / 2 . 

Therefore all the conditions are satisfied apart from the condition of continuity. /1 will be 
discontinuous when 

1 = ct(3xi + p) -1 / 3 , 3xi = a 3 — p. 
In order that this discontinuity coincides with the origin, it must be 

p = a 3 . (13) 

Therefore the condition of continuity relates in this way the two integration constants p and a. 
The complete solution of our problem reads now: 

where the auxiliary quantity 

i ?= (3x 1 + p) 1 / 3 = (r 3 + a 3 ) 1 / 3 
5 



has been introduced. 

When one introduces these values of the functions f in the expression (9) of the line element 
and goes back to the usual polar co-ordinates one gets the line element that forms the exact solution 
of Einstein's problem: 

ds 2 = (1 - a/R)dt 2 - - , - R 2 (d$ 2 + sin 2 M(fi 2 ), R = (r 3 + a 3 ) 1 / 3 . (14) 
1 — a/R 

The latter contains only the constant a that depends on the value of the mass at the origin. 

§5. The uniqueness of the solution resulted spontaneously through the present calculation. 
From what follows we can see that it would have been difficult to ascertain the uniqueness from 
an approximation procedure in the manner of Mr. Einstein. Without the continuity condition it 
would have resulted: 



( 3xi+p )-4/3 _ ( r 3 + ;9) -4/3 

h l-a(3 Xl +p)"V3 i_ a ( r 3 + p )-i/3- 



When a and p are small, the series expansion up to quantities of second order gives: 



h — —I 
r 4 



a , p 
1 + - - 4/3^ 



This expression, together with the corresponding expansions of fi, f'3, fi, satisfies up to the same 
accuracy all the conditions of the problem. Within this approximation the condition of continuity 
does not introduce anything new, since discontinuities occur spontaneously only in the origin. 
Then the two constants a and p appear to remain arbitrary, hence the problem would be physically 
undetermined. The exact solution teaches that in reality, by extending the approximations, the 
discontinuity does not occur at the origin, but at r = (a 3 — p) 1 ^ 3 , and that one must set just p = o? 
for the discontinuity to go in the origin. With the approximation in powers of a and p one should 
survey very closely the law of the coefficients in order to recognise the necessity of this link between 
a and p. 

§6. Finally, one has still to derive the motion of a point in the gravitational field, the geodesic 
line corresponding to the line element (14). From the three facts, that the line element is homoge- 
neous in the differentials and that its coefficients do not depend on t and on (fi, with the variation 
we get immediately three intermediate integrals. If one also restricts himself to the motion in the 
equatorial plane (1? = 90°, d-d = 0) these intermediate integrals read: 

R 2 ^. = const. = c, (16) 
ds 

(1 — a/R)^- = const. = 1 (determination of the time unit). (17) 



From here it follows 



^ +R 2 (l~ a/R) = =V [1 - h(l - a/R)} 
d(fi J c z 



6 



or with 1/R = x 



( dx \ 1 — h ha 9 o , . 

UJ = -^- + - 5 -*-* 2 W. (18) 



If one introduces the notations: c 2 /h = B, (1 — h)/h = 2A, this is identical to Mr. Einstein's 
equation (11), loc. cit. and gives the observed anomaly of the perihelion of Mercury. 

Actually Mr. Einstein's approximation for the orbit goes into the exact solution when one 
substitutes for r the quantity 



R=^ + a^=r{l + p j 



1/3 



Since a/r is nearly equal to twice the square of the velocity of the planet (with the velocity of light as 
unit), for Mercury the parenthesis differs from 1 only for quantities of the order 10 -12 . Therefore r is 
virtually identical to R and Mr. Einstein's approximation is adequate to the strongest requirements 
of the practice. 

Finally, the exact form of the third Kepler's law for circular orbits will be derived. Owing to 
(16) and (17), when one sets x = 1/R, for the angular velocity n = d<j)/dt it holds 

n = cx 2 (l — ax). 

For circular orbits both dx/dcfi and d 2 x/d(j) 2 must vanish. Due to (18) this gives: 



1 — h ha 2 s ^ a 2 

= — t, 1 — ttX — x + ax , = — 2x + oax . 

c z cr c z 



The elimination of h from these two equations yields 

a = 2c 2 x(l — ax) 2 . 

Hence it follows 



,2 



a o a a 



n = —x 



2 2R 3 2(r 3 + a 3 ) ' 

The deviation of this formula from the third Kepler's law is totally negligible down to the surface 
of the Sun. For an ideal mass point, however, it follows that the angular velocity does not, as with 
Newton's law, grow without limit when the radius of the orbit gets smaller and smaller, but it 
approaches a determined limit 



n 



aV2~' 



(For a point with the solar mass the limit frequency will be around 10 4 per second). This circum- 
stance could be of interest, if analogous laws would rule the molecular forces. 
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